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bJOl Abstract 

^ . 

We consider the spectral stability of solitary wave solutions (f){x)e '"* to the nonlinear Dirac equation in any 

dimension. This equation is well-known to theoretical physicists as the Soler model (or, in one dimension, the Gross- 

\f^ ' Neveu model), and attracted much attention for many years. We show that, generically, at the values of ui where the 

Vakhitov-Kolokolov stability criterion breaks down, a pair of real eigenvalues (one positive, one negative) appears 

from the origin, leading to the linear instability of corresponding solitary waves. 

As an auxiliary result, we state the Virial identities ("Pohozhaev theorem") for the nonlinear Dirac equation. 
We also show that ±2Loi are the eigenvalues of the nonlinear Dirac equation linearized at (^(a;)e~"^*, which are 
embedded into the essential spectrum as long as \uj\ > m/3. This result holds for the nonlinear Dirac equation with 
ri , any nonlinearity of the Soler form ("scalar-scalar interaction") and in any dimension. 

As an illustration of the spectral stability methods, we revisit Derrick's theorem and sketch the Vakhitov-Kolokolov 
stability criterion for the nonlinear Schrodinger equation. 
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^ ; 1 Introduction 

>■ 

fT^ ' Field equations with nonlinearities of local type are natural candidates for developing tools which are then used for the 

\C , analysis of systems of interacting equations. Equations with local nonlinearities have been appearing in the Quantum 

C^^ ' Field Theory starting perhaps since fifties ||Sch51a||Sch51bl . in the context of the classical nonlinear meson theory 

of nuclear forces. The nonlinear version of the Dirac equation is known as the Soler model IISol70l . The existence 
C^^ ' of standing waves in this model was proved in OSolTOI ICV86I . Existence of localized solutions to the Dirac-Maxwell 

^— ^ , system was addressed in ||Wak66l [Lis95l and finally was proved in IIEGS96I (for lu e (— m, 0)) and ||Abe98l (for 

u) e (— m, m)). The local well-posedness of the Dirac-Maxwell system was considered in llBou96| . The local and 
global well-posedness of the Dirac equation was further addressed in IIEV97I (semilinear Dirac equation in n = 3), 
BBouOOl (Dirac - Klein-Gordon system in n = 1), and in IIMNNO05I (nonlinear Dirac equation in n = 3). The 
\^ , question of stability of solitary wave solutions to the nonlinear Dirac equation attracted much attention for many 

^ ' years, but only partial numerical results were obtained; see e.g. IIAC81llAKV83llA"S83IIAS86llChu07l . The analysis 

of stability with respect to dilations is performed in ||SV86||CKMS10| . 

Understanding the linear stability is the first step in the study of stability properties of solitary waves. Absence 
of an eigenvalue with a positive real part will be referred to as the spectral stability, while its absence as the spectral 
(or linear) instability. After the spectrum of the linearized problem for the nonlinear Schrodinger equation IIVK73I 
was understood, the linearly unstable solitary waves can be proved to be ( "nonlinearly", or "dynamically") unstable 
IIGri88l IGOIOI . while the linearly stable solitary waves of the nonlinear Schrodinger and Klein-Gordon equations 
IISha83l [SS85I IWei86 1 and more general U(l)-invariant systems IIGSS87I were proved to be orbitally stable. The 
tools used to prove orbital stability break down for the Dirac equation since the corresponding energy functional is 
sign-indefinite. On the other hand, one can hope to use the dispersive estimates for the linearized equation to prove 
the asymptotic stability of the standing waves, similarly to how it is being done for the nonlinear Schrodinger equation 
IIWei85l , IISW92I , IIBP93I , ||SW99I , and QCucOll . The first results on asymptotic stabiUty for the nonlinear Dirac 
equation are already appearing flPSlOllBClll . with the assumptions on the spectrum of the Unearized equation playing 
a crucial role. 



C^ 



In this paper, we study the spectrum of the nonhnear Dirac equation hnearized at a sohtary wave, concentrating on 
bifurcation of real eigenvalues from A = 0. 

Derrick's theorem 

As a warm-up, let us consider the linear instability of stationary solutions to a nonlinear wave equation. 



■i/j = -Ai/j + giip), ip ^ Tpix^t) e 



X e 



n > 1. 



(1.1) 



We assume that the nonlinearity g{s) is smooth. Equation dl.lb is a Hamiltonian system, with the Hamiltonian 
m, ^) = /k. (4 + ^ + G(vi)) dx where G(s) = J^ gis') ds' . 

There is a well-known result IIDer64ll about non-existence of stable localized stationary solutions in dimension 
n > 3 (known as Derrick's Theorem). If M(a;, t) = 6{x) is a localized stationary solution to the Hamiltonian equations 
TT" = —5^E, ip = St^E, then, considering the family 6\{x) = 9{Xx), one has d\\^^-^E{(j)\) = 0, and then it follows 
that d1\^^-^E{(j)\) < as long as n > 3. That is, 6^E < for a variation corresponding to the uniform stretching, 
and the solution 6{x) is to be unstable. Let us modify Derrick's argument to show the linear instability of stationary 
solutions in any dimension. 

Lemma 1.1 (Derrick's theorem for n > 1). For any n > 1, a smooth finite energy stationary solution 6{x) to the 
nonlinear wave equation is linearly unstable. 



Proof. Since 6 satisfies -A9 + g{0) = 0, we also have -Ad^^O + g'{9)dx^9 



0. 



0, 



Due to lim 6{x) 

dxi9 vanishes somewhere. According to the minimum principle, there is a nowhere vanishing smooth function x S 
H°° (R") (due to A being elliptic) which corresponds to some smaller (hence negative) eigenvalue of L = —A + g'{9), 
Lx = ~(?X^ with c > 0. Taking ip{x,t) — 9{x) + r{x,t), we obtain the linearization at 9, —r = — Lr, which we 



rewrite as dt 


r 
s 


= 


■ 
-L 


1' 





r 
s 



The matrix in the right-hand side has eigenvectors 



X 

±cx 



corresponding to the 



eigenvalues ±c G M; thus, the solution 9 is linearly unstable. 

Let us also mention that d'^\^^„E{9 + tx) < 0, showing that S^E{9) is not positive-definite. 



D 



Remark 1.2. A more general result on the linear stability and (nonlinear) instability of stationary solutions to ( ILll i is 
in IIKS07II . In particular, it is shown there that the linearization at a stationary solution may be spectrally stable when 
this particular stationary solution is not from H^ (such examples exist in higher dimensions). 



Vakhitov-Kolokolov stability criterion for the nonlinear Schrodinger equation 

To get a hold of stable localized solutions. Derrick suggested that elementary particles might correspond to stable, 
localized solutions which are periodic in time, rather than time-independent. Let us consider how this works for the 
(generalized) nonlinear Schrodinger equation in one dimension. 



^^ti: = -\^li: + gM'')^l^, 



ip = ip{x,t) e C, x G 



t e 



(1.2) 



where g{s) is a smooth function with m := g(0) > 0. One can easily construct solitary wave solutions (j){x)e^^'^* , 
for some w G M and (f) G H^{M.): (f){x) satisfies the stationary equation ujcj) = ^\4>" + 9{4>'^)4>^ ^nd can be chosen 
strictly positive, even, and monotonically decaying away from x = {). The value of u can not exceed m. We consider 
the Ansatz il^i^i t) = {4'{x) + p{^i t))e~'^'^*, with p{x, t) G C. The linearized equation on p is called the linearization 
at a solitary wave: 

Bcp{x,t) 



5tR = JLR, 



R{x,t) 



Jnip{x, t) 



with 



1 
-1 



L_ 



L_ = 



2 -^ 



L+ = L_ + 2g'{ 



(1.3) 



(1.4) 



Note that since I 7^ 1-+, the action of L on p considered as taking values in ' 

lim 0(x) = 0, the essential spectrum of I and L_|- is [m — lj, +cxd). 



is R-linear but not C-linear. Since 



\x\ 



First, let us note that the spectrum of JL is located on the real and imaginary axes only: o'(JL) C 



Uil 



To 



prove this, we consider (JL)^ 











Since I is positive-definite (cj) G ker I , being nowhere zero. 



corresponds to its smallest eigenvalue), we can define the selfadjoint root of I ; then 

a,((JL)2)\{0} = a,(L_L+)\{0} = a,(L+L_)\{0} = a.iL'J'UL'J'UO} C M, 

with the inclusion due to L_ L-|_L_ being selfadjoint. Thus, any eigenvalue A e (^^(JL) satisfies A'^ e R. 

Given the family of solitary waves, 0^(a;)e"*"*, w G fi C K, we would like to know at which ui the eigenvalues 
of the linearized equation with Re A > appear Since A^ G K, such eigenvalues can only be located on the real axis, 
having bifurcated from A = 0. One can check that A = belongs to the discrete spectrum of JL, with 



JL 







0, 



JL 



-d, 







for all uj which correspond to solitary waves. Thus, if we will restrict our attention to functions which are even in 
X, the dimension of the generalized null space of JL is at least two. Hence, the bifurcation follows the jump in the 
dimension of the generalized null space of JL. Such a jump happens at a particular value of lj if one can solve the 



equation JLa 







This leads to the condition that 







is orthogonal to the null space of the adjoint to 



JL, which contains the vector 







; this results in {(j)^, duj4>uj) = duj\\4>uj\\i2/'2 = 0. A slightly more careful analysis 
IICP03I based on construction of the moving frame in the generalized eigenspace of A = shows that there are two 



real eigenvalues ±A G M that have emerged from A = when uj is such that 9cj||'/'aj||?2 becomes positive, leading to 



a linear instabiUty of the corresponding solitary wave. The opposite condition. 



d^ 



"t^WL^ 



<o, 



(1.5) 



is the Vakhitov-Kolokolov stability criterion which guarantees the absence of nonzero real eigenvalues for the nonlin- 
ear Schrodinger equation. It appeared in IIVK73I ISha83 GSS87I in relation to linear and orbital stability of solitary 
waves. The above approach fails for the nonlinear Dirac equation since L_ is no longer positive-definite. 

For the completeness, let us present a more precise form of the Vakhitov-Kolokolov stability criterion IIVK73I . 

Lemma 1.3 (Vakhitov-Kolokolov stability criterion). There is A G (Tp(JL), A > 0, where JL is the linearization (11.31 1 

at the solitary wave (p^{x)e~'^'^^, ifandonly if -j^\\4'ui\\'j;^2 > at this value ofuj. 



Proof. We follow IIVK731 . Assume that there is A G crd(JL), A > 0. The relation (JL - A)S = impUes that 
A^Si = —I L-|_Si. It follows that Si is orthogonal to the kernel of the selfadjoint operator I (which is spanned by 

i,Si) = -4(<^,-L-L+Si) = -^(L_0,-L+Si) = 0, 



hence there is ?] G L^ 



1— 



A2 
^, C) such that Si 



A^ 



I f] and X^T] 



-L+Si. Thus, the inverse to I can be applied: 



A^LI^Si 



-L_|_Si. Then 



A2(ry,L_r;) = -(Si,L+Si). 



Since I is positive-definite and 77 ^ ker I , it follows that {?], I rj) > 0. Since A > 0, (Si, L+Si) < 0, therefore 

the quadratic form (■, L+-) is not positive-definite on vectors orthogonal to 0^^. According to Lagrange's principle, the 
function r corresponding to the minimum of (r, L+j-) under conditions (r, (f>u:) ~ and (r, r) = 1 satisfies 



Ur = 



a, 13 e 



(1.6) 



Since (r, L_|_7-) = a, we need to know whether a could be negative. Since L^dx<pui = 0, one has Ai = G (Tp(L_|_). 
Due to dx(t>ui vanishing at one point (x = 0), there is exactly one negative eigenvalue of L+, which we denote by 



•^0 £ cTp (!-+)• (This eigenvalue corresponds to some non-vanishing eigenfunction.) Note that /3 7^ 0, or else a would 
have to be equal to Ao, with r the corresponding eigenfunction of L+, but then r, having to be nonzero, could not be 
orthogonal to 0CJ . Denote A2 = inf((T(L+) nR+) > 0. Let us consider /(z) — {(j)uj,(L+ — z)~^(f>i^), which is defined 
and is smooth for z G (Ao, A2). (Note that f{z) is defined for z = Ai = since the corresponding eigenfunction dx(l>uj 
is odd while (f)u: is even.) If a < 0, then, by ( 11.6b . we would have /(a) ~ {(f)u:, (L+ — a)^^(f)ui) = jji't'uj, '') = 0, and 
since f'{z) > 0, one has /(O) > 0. On the other hand, /(O) = (<?i<^, L^V^) = {(f>u:,d^(l)^) = i^ /^ |0^(a;)P dx. 
Therefore, the linear instability leads to a < 0, which results in -^ J^ \4'lj{x)\'^ dx > 0. 

Alternatively, let ^||(?!'w|li2 > 0. We consider the function /(z) = (^(^,(1+ — z)~^(j}i^), z e p(l-+). Since 
/(O) = ((/)„, L^^cj)^) > 0, /'(z) > 0, and lim /(z) = —00 (where Aq < is the smallest eigenvalue of 1+), there 

2-i-Ao+' 

is a 6 (Ao,0) C /o(L+) such that /(a) = {(j)^, (L+ — a)^^(f)^) = 0. Then we define r = (L+ — a)^^(f)ui. Since 

{(pi^, r) = f{a) ~ 0, there is 77 such that r = I 77. It follows that the quadratic form L_ 1-+L_ is not positive 

definite: 



(L^?7, (L^L+L^)L: 
Thus, there is A > such that -A^ e a{\}!'^L+\}!'^)\ then also -A^ £ (t(L_L 



eigenvector, I L_|-^ = — A^^; then 



f]) = {r, l+r) 
1/2 







iwe 



(r, {ar + (j)^)) = Q:(r, r) < 0. 

Let ^ be the corresponding 
hence A G cr(JL). D 



= A 






^-L+e 



Our conclusions: 

1. Point eigenvalues of the linearized Dirac equation may bifurcate (as w changes) from the origin, when the 
dimension of the generaUzed null space jumps up (when the Vakhitov-Kolokolov criterion breaks down). 

2. Since the spectrumof the linearization does not have to be asubsetof KUiM, there may also be point eigenvalues 
which bifurcate from the imaginary axis into the complex plane. (We do not know particular examples of such 
behavior for the nonlinear Dirac equation.) 

3. Moreover, there may be point eigenvalues already present in the spectra of linearizations at arbitrarily small 
solitary waves. Formally, we could say that these eigenvalues bifurcate from the essential spectrum of the free 
Dirac operator (divided by i), which can be considered as the linearization of the nonlinear Dirac equation at the 
zero solitary wave. 



In the present paper we investigate the first scenario. The main result (Lemma 14.1b states that if the Vakhitov- 
Kolokolov breaks down at some point a;*, then, generically, the solitary waves with lj from an open one-sided neigh- 
borhood of oj* are linearly unstable. 

We also demonstrate the presence of the eigenvalue ±.2uji in the spectrum of the linearized operator (Corollarv l2.8b 
and obtain Virial identities, or Pohozhaev theorem, for the nonlinear Dirac equation (Lemma [3.2b . which we need for 
the analysis of the zero eigenvalue of the linearized operator. 



2 Linearization of the nonlinear Dirac equation 

The nonlinear Dirac equation in M" has the form 






(2.1) 



where dj = -^, N is even and g smooth, with m := g(0) > 0. The Dirac matrices aj and /3 satisfy the relations 

ajak + ataj = 2SjklN, aj(3 + (3aj ^0, (3'^^ In, 'i^ < j,k <n, 



where In is an N x N unit matrix. We will always assume that f3 

cr 



N/2 









-/, 



Af/2_ 



In the case n ~ 1, we assume 



«! = —0-2; in the case n = 3, one could take a, 



a, 



, where aj are the standard Pauli matrices. Equation 



( 12.1b . usually with g{s) = 1 — s, is called the Soler model ISolTOI , which has been receiving a lot of attention in 
theoretical physics in relation to classical models of elementary particles. 

Remark 2.1. In terms of the Dirac 7-matrices, equation (12. Il l takes the explicitly relativistically-invariantform 17'' 5^^ = 
g{'ip*/3ip)ip, where 7" = (3, 7^' = Pa,, do = dt, dj = d^^. 



Definition 2.2. The solitary waves are solutions to ( 12. Il l of the form 0(^(a;)e *"*, ^^^ G H^il 
Below, we assume that there are solitary waves for u) from some nonempty set 17 C R: 



,C 



N\ 



iO £ 



0^6 



w e r2 c 



eil^( 



^N 



), 



(2.2) 



with (f),^ smoothly depending on to. 

We will not indicate the dependence on uj explicitly, and will write 4> instead of 4>^. 
The profile of a stationary wave satisfies the stationary nonlinear Dirac equation 

n 

C-(j) := {-i^ajdj - uj + g{4>* P(I))P)(I) = 0. 
i=i 

The energy and charge functionals corresponding to the nonlinear Dirac equation (12. Il l are given by 



(2.3) 



/^„ I -»EV'*«j5,V; + G(V'*/3^) rx, 



where G{s) is the antiderivative of g{s) which satisfies G{0) = 0. Q{ip) is the charge functional which is (formally) 
conserved for solutions to ( 12. Il l due to the U(l)-invariance. The nonlinear Dirac equation ( 12. Il l can be written in the 
Hamiltonian form as 9t Im?/; = —^Sbc^E, dt^ip = ^(5inn/;£', or simply ip ~ —i5.^'E. The relation (12.31 ) satisfied 
by the profile of the solitary wave 0(a;)e~'"* can be written as 



E'{4>) = ujQ'{4>), 



(2.4) 



where the primes denote the Frechet derivative of the functionals E{ip), Q{ip) with respect to (Re ^, Ln ?/j). 

Let us write the solution in the form ip{x, t) — {(j>{x) + p{x, i))e~*'^*, p{x, t) G C^. The linearized equation on p 
is given by 

P = JCp, (2.5) 

where J corresponds to a multiplication by 1/i and 

Cp = C-p + 2g'{4)*l34))l34)Y{e{(t)*Pp). 

Note that, because of the presence of 'Rc{(t>* jip), the action of £ on p G C^ is M-linear but not C-linear. Because of 
this, it is convenient to write it as an operator L acting on vectors from R^^; then (I2.5l l takes the following form: 





9iR 


= JLR, L = 


JA,5, 


where g = 


= ff(r/3</'), 9' = .9'(r /30) and 




<? = 


■Re0 
Ltk/) 


, J- 




-In 



In 




,gB + 2g'B<Z'(B<?, •)R2iv, 



Re a j — Im a^ 
Ima, Reai 



R(a;,i) 



B 



Rep 
Imp 



'(i 0' 
/3 



(2.6) 



(2.7) 



Note that J, Aj, and B correspond to multiplication by — i, aj, and /? under the C ^-> M correspondence. 
Remark 2.3. When n = 1, one can take ai = — (T2 (so that A^ = 2). Then L has a particularly simple form since 
G C^ can be chosen valued in R^: L= + , with L_ = icra^x - tj +5/3, L+ = L_ + 2g'/30 (/30, ■ )c2. The 

numerical and analytical study of spectra of I , L-|_ in this case is contained in IIBC09I . 



Lemma 2.4. o-ess(L) = Uessi^) = R\(— m — a;,m — w); (Tess(JL) = i(K\(— m + tj,m — w)). 

Proof. One has a{—iajdj + l3m) = K\(— m, m). Note also that {—iajdj + j3mY = —A + m? has a spectrum 
[m? , oo). Taking into account that the symbol of C at |a;| -^ c« is aj^j + (3in — uj, one concludes that (Tess(i-) = 
(Jessi^) = M\(— m ~ uj,m — uj). Since the eigenvalues of J are ±i, corresponding to clock- and counterclockwise 
rotations in C, one deduces that cTess (JL) = i(K\(— m + lj, m — a;)). D 

Lemma 2.5. The null space of Jl. is given by ker JL = Span {J<P, dk'P] 1 < k < n\ . 

Proof Recall that 

n n 

Since 0(a;) e C^ satisfies the stationary nonlinear Dirac equation (|2.3| i. we get: 

/:(-i0) == C- {-i(j)) + 2 Rc((/)*/3(-i0)) = 0. (2.8) 

Taking the derivative of ( 12.31 ) with respect to Xk yields 

n 

D 

Lemma 2.6. Let a^ be an hermitian matrix anticommuting with aj, 1 < j < n, and with /?. Then aQ(j) is an 
eigenfunction of C- and of C, corresponding to the eigenvalue A ~ —2u. 

Remark 2.7. If n = 3, one can take ao = Qia2a3/3. 

Proof. Since ag anticommutes with aj (1 < j < n) and with /3, and taking into account (12.3b . we have: 

n n 

£-.ao(f> — {~i y^ ctjdj — oj + g{(f>* l3(f))[3)ao(f> — ao{i /^ ctjdj —uj — g{(f>*l3(f>)l3)(j) — ao(— £- -~2uj)(f) = —2ujao(j). 



Since ao and (3 are Hermitian, 2Re[(/)*/3ao0] = (p* l3ao(f> + (j)*f5aQ(j) = 0*{/3,q;o}0 ~ 0; therefore, one also has 
Cao(t) — C^aQcf) = —2u)a()(f). D 

It follows that the linearization operator has an eigenvalue 2ujv. 

2uji G (T.p{JC) = <Tp(JL). 

Since fT(JL) is symmetric with respect to M and iR, for any g{s) in (12.1b and in any dimension n > 1, we have: 

Corollary 2.8. ±2a;i are L^ eigenvalues of]\-. 

Remark 2.9. For \uj\ > m/i, the eigenvalues ±2uji are embedded in the essential spectrum. This is in contradiction 
with the Hypothesis (H:6) in IBCllI on the absence of eigenvalues embedded in the essential spectrum, although we 
hope that this difficulty could be dealt with using a minor change in the proof. 

Remark 2.10. The result of Corollarv l2.8l takes place for any nonlinearity g{ip* (3ip) and in any dimension. The spatial 
dimension n and the number of components of ijj could be such that there is no matrix ao which anticommutes with aj , 
1 < j < n, and with /3; then the eigenvector corresponding to ±.2uji can be constructed using the spatial reflections. 



3 Virial identities 

When studying the bifurcation of eigenvalues from A = 0, we will need some conclusions about the generalized null 
space of the linearized operator. We will draw these conclusions from the Virial identities, which are also known as the 
Pohozhaev theorem ||Poh65l . In the context of the nonlinear Dirac equations, similar results were presented in IIES95I . 

Lemma 3.1. For a dijferentiable family \^^ 4)x E H^{W^), (f)\\^^^ ^ (p, one has d\E{(j>\)\^^-^ = ujd\Q{(j)\)\^^-^. 

Proof. This immediately follows from jlAj . D 

n 

We split the Hamiltonian E{4>) into E{^) = T{4>) + V{^jj) = J2 TjW + V{^p), where 

i=i 

T,-(^) = -i/" tl;*ajdji!(Px (no summation in j), T{iP) ^y^Tj{i,), ¥{%{;) ^ [ G{f* (ii^)(r'x. 

Lemma 3.2 (Pohozhaev Theorem for the nonlinear Dirac equation). For each solitary wave (^(a;)e~*"*, there are the 
following relations: 

n ~ 1 If 

r(^) + 1/(0) = wQ(0), T,{^)^-T{cf>)^ {G{4>*p4>)~<l^*P<t>9{<l^*P<ty))dx, l<j<n. 



Proof We set 0a (x) = 0(.ti/A, .T2, . . . , a;„). Since Ti(</.a) = Ti(0), r,(0A) = AT,(0) for 2 < j < n, V{(jyx) = 
Ay (0), Q{4>\) = \Q{4)), we can use Lemma [TTI to obtain the following relation ("Virial theorem"): 

n n 

^T,(0)+F(0)=9aL.,(E^:>(<^a) + ^('^a)) =c.9AL.,g(0A)-c.Q(0). (3.1) 

Similarly, rescaling in Xj, 1 < j < n, we conclude that T'i(0) = . . . = Tn{(j)) = —T{(j)), and ( 13.1b gives 

n — 1 



'-T{(f) + V{(f)^^Q{4,). (3.2) 

li 

Moreover, the relation ( 12.3b vields T((/)) + /g„ (j)* jicj) g[cj)* jicj)) dx = ujQ{(f>). Together with ( 13.21 ). this gives the desired 

relation iT(0) = 4„ (G(0*/30) - 0*/305(0*/30)) dx. D 

Remark 3.3. For all nonlinearities for which the existence of solitary wave solutions is proved in IIES95II . one has 
G{s) — sG'{s) > for all s £ K except finitely many points (e.g. s = 0); hence for these sohtary waves one has 
T(0) > 0. In particular, for G{s) :^ s- ^,one has T{(j)) ^ n /g„ I^IM! dx > 0. 

4 Bifurcations from A = 

Lemma 4.1. Assume that the nonlinearity satisfies the following inequality (see Remark \3.3\l : 

77 + 1 

——G{s)~sG'{s)>0, seR, s^O. (4.1) 

n 

Further, assume that <jf (p and (f)* P4> are spherically symmetric and that 

.yV{V} = Span{ J<?, 9fc<?; 1 < fc < n} ; dim^(L) = n+ 1. 

IfdujQ{4') ¥" 0> f/^^" f^^ generalized null space of]l. is given by 

^(JL) = Span { ]<P, d^^, dk^, Ak<P\ l<k< n} ; dim^(JL) = 2n + 2. 

IfdujQ{4i) vanishes at a;*, then dim,/(^(JL|^ ) > 2n + 4. Moreover, generically, there is an eigenvalue A G crrf(JL) 
with Re A > for lo from an open one-sided neighborhood of ui^,. 



Remark 4.2. The assumption that 0*c/) and (t>* Pi' ^re spherically symmetric is satisfied by the ansatz 



cj,{x) 



1 



cosf 



used in e.g. IIWak66l ISol70l IES95I IEGS961 . 

Proof. Taking the derivative of ( I2.3l l with respect to oj, we get 

n 



codu^d) - = -^9,^0 -0 = 0. (4.2) 



Since 4>* f3ak4> + {oik4>)* P<i> = 0*{/55Q;}0 = 0, we have 

C{ak4') = C^{ak(i)) = -2idk<p - Swa^^ - a^C^cp = -2idk<p - 2ujak4i- (4-3) 

Similarly, since (p* (3{ixk<p) + {ixk4>)* P(t> — 0, 

C{ixk4>) = C- {ixk(t>) = ak(t> + ixkC-(t> = afe0. (4.4) 

Using ( |4.31 l and (|4.4| i. we have 

C{ak<i) + 2ujixk<i)) = -20fe0. (4.5) 

Until the end of this section, it will be more convenient for us to work in terms of ^ G K^^ (see ( I2.7l i). We summarize 
the above relations dZST l, (|Z9] |. (|42]i, ( l43T l. (|4!4] |, and ( l43T l as follows: 

LJ<5 = 0, Lafc<? = 0, Ld^$^$, U-Xk^^Ak^, UAk^-2ioxkJ$)^2Jdk^, 



L = JAj9j - tj + gB + 2g'B<P{<P, B • )k2« , Lq = JAja^- - w + gB, 



where 



with g = g{(t)*li(t)) = gi^*'B<P), g' = g'{<f)* Pep) = g'(<?*B<?). There are no Ffc such that JLFfc = Afc<P - 2ujXk]<P. 
Indeed, checking the orthogonality of A^^ — 2ujXk]'P with respect to ,yK((JL)*) = Span{^, J9fc<?}, we have: 



{Ak^~2uxkJ$,Jdk<P)^ f 4>*{^iy2^kdu)4>dx + Lo j (b*(bdx 



r(0) 



0.0(0). 



Note that there is no summation in k in ( 14.6b . By Lemma [3^ the right-hand side of ( I4.6l l is equal to 

- + cjQ = - + l!l:illl + T/ = T + V^= / (n(G(p) - pG'(p)) + G{p)) dx, 
n n n Ju- 

where p{x) — (p* (x)/30(x). By (14. It , the right-hand side of ( 14. 7b is strictly positive. 



(4.6) 



(4.7) 



Due to Lemma 12.41 we can choose a small counterclockwise-oriented circle 7 centered at A = such that at 
uj = uj^, the only part of the spectrum cr (JL) inside 7 is the eigenvalue A = 0. Assume that O is an open neighborhood 
of w* small enough so that (t(JL) does not intersect 7 for u e O. Define 



Po 



dX 



27riJ^X-jV 



wgC 



(4.8) 



For each ui G O, Po is a projection onto a finite-dimensional vector space ^ := Range(Po) C i^(M"). The 
operator PqJI- is bounded (since Pq is smoothing of order one) and with the finite-dimensional range. Applying the 
Fredholm alternative to P0JI-, we conclude that there is E^ such that P0JI--E3 = d^j'P (hence JL£^3 = dui'P) if and 



only if dujQ{(j>) ~ 0. Indeed, one can check that it is precisely in this case that dui'P is orthogonal to c/K((JL)*) 
Sp&n{^Jdk^}: 

The right-hand side vanishes since it is the fcth component of the (zero) momentum of the standing solitary wave. 
Remark 4.3. Using ( 12. 3t , one can explicitly compute 



2uj{4>,ak4>) = {4>,ak{-iajdj + gf^)<l)) + {{-iajdj + gl3)(j),ak(l)) = -2i{4>,dk4>) = ^M 



dM*(IAdx^O. 



Once there is £^3 such that JLi^a = d^^P, there is also E^ such that JL£J4 = £^3, since E^ is orthogonal to the 
null space ^((JL)*) = Span{<?, jafc<?}: 

(£;3,<?) - {E:i,w^^) = (L£;3,a^<P) = -{]d^$,d^$) = o, 

2(£;3, J5fc<?) = (£;3,L(Afc^ - 2uxu]^)) = (L£;3,Afc<? - 2c^.TfeJ<?) = -(ja^<P, Afe-P - 2wxfcJ<Z') = 0. 
To check that the right-hand side of the second line is indeed equal to zero, one needs to take into account the following: 

(J5^<?, ujXk']'P) - oj{d^'P, Xk'P) - ^9^ / xu(l)*(l) dx = 0, (4.9) 

{]d^$,Ak<P) = (ja^<l',L(-a-J<?)) = {]d^^,U{-XkW)) = -{Ud^^,Xk^) 

= ~{<P ~ 2g'icl)* l3(l>)B<P{<P,Bd^^)^2N ,Xk'P) = -{'P,Xk'P) + {2g'{(j,* (3(l))B<P{'P,Bd^'P)^2N ,Xk<P) 

Xk(j)*(l)dx + d^ f XkK{(j)*P(f))dx^Q. (4.10) 

Above, K{s) is the antiderivative of sg'{s) such that /v (0) — 0. The integrals in the right-hand sides of ( 14.91 ) and 
(14.10b are equal to zero due to our assumption on the symmetry properties of cf)* (f) and 0* f3(j). Note that £_ is C-linear, 
hence commutes with a multiplication by i, and therefore J and Lo commute; we used this when deriving ( 14.101 ). 
We will assume that 

(£;3,L£;3)7^o. (4.11) 

Then there is no £^5 such that ]\-E^ ~ E4, since E4 is not orthogonal to the null space J\f{(J'L)*) 3 <P: 

{Ei,$) = {Ei,\.d^$) = (£4,LJL£3) = -{E^,\-E:i) + 0. (4.12) 

Remark 4.4. If (14.11b is not satisfied, then the dimension of the generalized null space of JL at cj* may jump by more 
than two; this means that there are more than two eigenvalues colliding at A = as w passes through w*. We expect 
that generically this scenario does not take place. 

We will break the finite-dimensional vector space ^ = Range (Pq) into a direct sum 

SC ^'3^ ®^, S" = Span {dj<l>, Aj-P; I < j < n} , 

so that both ^ and 2f are invariant with respect to the action of JL. (Let us mention that Pq, ^, '3/ , ^, and JL 
depend on w.) Set 

64(0;) = P0-E4 e ^, 63(0;) = JLe4 e ^. (4.13) 

Since Pq continuously depends on lo, 63(0;) and 64(0;) are continuous functions of cj. {e\{uj\ 62(0;), 63(0;), 64(0;)}, 
cj G O, is a frame in the space '3^ . For some continuous functions cri(cj), 0-2 (cj), a-j,(LS), and a^^iuj), there is the 
relation 

JLe3(u;) = o-i(a;)ei(a;) + o-2(a;)e2(a;) + cr3(u;)e3(a;) + cr4(u;)e4(a;), w e O. (4.14) 



Evaluating (I4.14l i at w* , we conclude that 

cri(a;*)=0, a2{u}*) = I, cr3(w*) = 0, (74(w*)=0. (4.15) 

In the frame {ei(tL)), 62(0;), 63(0;), 64(0;)}, the operator JL restricted onto ^ is represented by the matrix 



M 






1 


CTl(w) 











CT2(w) 











CT3(CJ) 


1 








CT4(w) 






ujeO, 



(4.16) 



which is a 4 X 4 Jordan block at w*. Let us investigate its entries. Pairing ( |4.14| i with <P = J^^ei and taking into 
account that (j"^ei, (JL)2e4(w)) = - (Lei, JLe4(w)) = 0, (J"^ei,ei) = 0, and (J^'^ei,e3) = (J^^ei, JLe4) = 
— (Lei, 64} = 0, we get: 

= a^iio) (j-'ei, ea) + ai{Lo) (j-^ei, e4) . (4.17) 



Define the function 
Since 



^(w) :=-(<?, e4). 



ujeO. 



fJ,{uJ^, 



-(<?,e4)L^ =-(<?L^,£;4) = (£;3,L£;3)7^o 



(4.18) 



(4.19) 



by ( |4.12t . we may take the open neighborhood O of w* to be sufficiently small so that fJ,{uj) does not vanish for w e C 
Then 

0-4 (w) = ^^; , LoeO. (4.20) 

Ai(cj) 

Let us show that 173(0;) is identically zero in an open neighborhood of oj* . Applying (JL)'^ to ( 14.14b . we get: 



(JL)3e3(a.) = a3ioj)OLfe3iu) + a4{Lo)OLfei{uj). 

Coupling with J^^e4 and using (J^^e4, (JL)^e3) = (e3,LJLe3) = and (J^^e4, (JL)^e4) 
(due to anti-selfadjointness of LJL), we have 



(4.21) 
(e4,LJLe4) =0 



a3(w)(rie4,(JL)2e3)-0. 
The factor at 173(0;) is nonzero. Indeed, using (14.151 1. 

(J"^e4,(JL)2e3)L^ = (J"ie4, cr2ei+CT3jLe3 + cr4e3)L^ =(rie4,ei)L^ =-(e4,<?)L^, 



(4.22) 



which is nonzero due to (|4.12t . We conclude that <t-^{uj) is identically zero in an open neighborhood of w* . We take O 
small enough so that 0-3 j^ =0. 

Near u = ui^, the eigenvalues of JL which are located inside a small contour around A = coincide with the 
eigenvalues of the matrix M, defined in ( 14.161 1. Since 0-3 is identically zero in O, these eigenvalues satisfy 



X^X^ - aiiiu)) = 0, 



weO. 



(4.23) 



By (I4.2OI 1. if di^Q{4)) changes sign at a;*, so does 174(0;), hence in a one-sided open neighborhood of o;* there are two 
real eigenvalues of JL, one positive (indicating the linear instability) and one negative. D 



Remark 4.5. This argument is slightly longer than a similar computation in IICP03I since we do not assume that (p could 
be chosen purely real (allowing for a common ansatz used in IIES95I in the context of the nonlinear Dirac equation), 
and consequently we could not take Bj to be "imaginary" for j odd and "real" for j even, and enjoy the vanishing of 

(J^^ej, efc) for j + fc is even. 
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5 Concluding remarks 

In the conclusion, let us make several observations. 

Remark 5.1. In the case of the nonlinear Schrodinger equation, the function yu(a;) is strictly positive (this is due to 

positive-definiteness of 1 in ( 11.41 ). which results in positivity of ( |4.19t ; see IIVK73I ICP03I ). so that the collisions 

of eigenvalues at A = always happen according to the following scenario: when dQ/du! changes from negative to 
positive (no matter whether this happens as uj increases or decreases), there is a pair of eigenvalues on the imaginary 
axis colliding at A = and proceeding along the real axis. Then, further, if dQ/duj changes from positive to negative, 
this pair of real eigenvalues return to A = and then retreat onto the imaginary axis. For the Dirac equation, we 
can not rule out that /i(a;) changes the sign (becoming negative). If this were the case, vanishing of dQ/dtu would 
be accompanied with the reversed bifurcation mechanism: as dQjduj goes from positive to negative, another pair 
of imaginary eigenvalues collide at A = and proceed along the real axis. We do not have examples of particular 
nonlinearities which lead to such a scenario. 

Remark 5.2. In one dimension, for the nonlinearity g(.s) = 1 — s*^', fc G N, using the asymptotics as in IIGuaOSI . one 
can derive that, as cj ^> 1—, the function /i(w) defined in ( |4.18t has the asymptotics 

Ai(w) = 0((l-cj2)i:-5). 

There is the same asymptotics in the case of the nonlinear Schrodinger equation (11.21) with the same nonlinearity 
g{s). In particular, lim /i(w) = +oo, hence /i(w) remains positive for lo sufficiently close to 1 (precisely as for the 

u;— fl — 

nonlinear Schrodinger equation, when /^t(cj) > for all uS), suggesting that the bifurcation scenario for the nonlinear 
Dirac equation near w = 1 is the same as for the nonlinear Schrodinger equation. 

Remark 5.3. For the nonlinear Dirac equation in one dimension with the nonlinearity gi^s) = 1 — s ("Soler model"), 
the solitary waves exist for uj G (0, 1). One can explicitly compute that the charge is given by Q(uj) = /^ |V'|^ Ax = 

^^"'^"'^ (see e.g. IILG75II ) so that dQ/dio < for uj e (0, 1), implying that there are no eigenvalues of JL colliding at 
A = for any uj € (0, 1). According to our numerical results IIBC09L we expect that in this model the solitary waves 
are spectrally stable, at least for uj sufficiently close to 1. 

Remark 5.4. Even if dQ/duj never vanishes, so that there are no bifurcations of nonzero real eigenvalues from A = 0, 
there may be nonzero real eigenvalues present in the spectrum of all solitary waves. We expect that the Vakhitov- 
Kolokolov criterion could again be useful here when applied in the nonrelativistic limit uj —i' m, when the properties 
of the nonlinear Dirac equation are similar to properties of the nonlinear Schrodinger equation. This idea has been 
mentioned in MCKMSIOI . Our preliminary results indicate that in one dimension, for uj sufficiently close to 1, the 
nonlinearity with g{s) = 1 — s'^ + o{s'') with k = 1 and fc = 2 does not produce nonzero real eigenvalues, while for 
fc > 3 there are two real eigenvalues, one positive (leading to Unear instability) and one negative. 
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